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The positron-neutrino correlation in the + — » + (3 decay of 32 Ar was measured 
at ISOLDE by analyzing the effect of lepton recoil on the shape of the narrow proton 
group following the superallowed decay. Our result is consistent with the Standard 
Model prediction; for vanishing Fierz interference we find a = 0.9989 ± 0.0052 ± 
0.0036. Our result leads to improved constraints on scalar weak interactions. The 
positron-neutrino correlation in 33 Ar decay was measured in the same experiment; 
for vanishing Fierz interference we find a = 0.944 ± 0.002 ± 0.003. The 32 Ar and 
33 Ar correlations, in combination with precision measurements of the half-lives, 
superallowed branching ratios and beta endpoint energies, will determine the isospin 
impurities of the superallowed transitions. These will provide useful tests of isospin- 
violation corrections used in deducing |Vud which currently indicates non-unitarity 
of the KM matrix. 
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1. The (e, u) Correlation as a Probe of Physics Beyond the Standard 
Model 

In the Standard Model, nuclear f3 decay is mediated by the exchange of W 
bosons which have only vector and axial-vector couplings. However, extensions 
of the standard model, such as super-symmetric theories with more than one 
charged Higgs doublet, or leptoquarks, naturally predict scalar or tensor weak 
couplings JTJ. A general effective Hamiltonian for allowed (3 transitions that 
respects Lorentz invariance is [0] 

+ (lpplpn)(Cs1pe1pu + C' s 1p e "/ 5 1p v ) 

+ Hermitian conj . (1) 

where a term proportional to (V'pTsV'n) has been neglected because nucleons 
are non-relativistic. In the Standard Model, Cy = C' v , Ca = C' A , and 
Cs = C' s = Ct = C' T = 0. Jackson et al. Q computed the nuclear-/?-decay 
rate from this Hamiltonian; for an unoriented initial state and summation over 
the lepton helicities 

dW = dW (l + + , (2) 

where 

at = \M F \ 2 (\C V \ 2 + \C' V \ 2 - \C S \ 2 - \C' S \ 2 ) 

-i|M GT | 2 (|^| 2 + \C' A \ 2 - \C T ? - \C' T \ 2 ) (3) 

b£ = ±2 7 Re[\M F \ 2 (C v C s + C'^C's) 

+ \M GT \ 2 (C A C T + C'XC' T )] (4) 

and 

£ = \M F \ 2 {\C V \ 2 + \C' V \ 2 + \C S \ 2 + \C' S \ 2 ) 

+ \M GT \ 2 (\C A \ 2 + \C' A \ 2 + \C T \ 2 + \C' T \ 2 ) , (5) 

where Mp and Mqt are rank-0 and rank-1 nuclear matrix elements. We have 
simplified the expressions by neglecting the small Coulomb effects, and adopting 
the allowed approximation (ignoring curvature in the lepton wave functions and 
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the recoil-order corrections). In what follows we shall assume Cy = C' v and 
express the scalar couplings in terms of: 

Cs = Cs/C v , and C' S = C' S /C V . (6) 

Precise measurements of the e-v correlation coefficient, a, and of the Fierz- 
interference term, b, can potentially yield information about physics beyond the 
standard model. Equations 2-5 are complicated and depend on nuclear physics 
information through Mp and Mqt- However, for pure Fermi or pure GT transi- 
tions the expressions simplify and the Mp and Mqt factors cancel. In pure Fermi 
transitions, a positron- neutrino correlation coefficient a < 1 would immediately 
imply the presence of scalar currents; while in GT decays a value a > —1/3 would 
imply tensor currents. 

To measure the e-v correlation one must determine the daughter's velocity, 
because it is out of the question to detect the low-energy neutrino. But, in gen- 
eral, measuring the daughter's recoil velocity is not a trivial task: the velocities 
are very small and the interaction of the daughter with the surrounding medium 
can jeopardize the measurement. Fortunately nature has provided us with cases 
where the daughter states are unbound to proton emission, so that the daughter's 
velocity can be determined via the 'Doppler' broadening of the beta-delayed pro- 
ton groups. Protons are preferred over other /3-delayed radiations because they 
are emitted before the daughter nucleus has slowed appreciably and can be de- 
tected with high resolution and good efficiency. Neutrons, on the other hand, are 
hard to detect with high resolution and high efficiency simultaneously. Gamma 
rays are generally emitted too slowly to probe the daughter velocity before it has 
been reduced by interaction with the surrounding medium. This paper discusses 
superallowed transitions from proton-rich nuclei with proton-unbound daughter 
states. These transitions are strong, which enhances the signal-to-noise ratio, and 
the widths of daughter states are ideally suited for probing the e-v correlation: 

• the time scale for proton decay is short enough so that the decays occur be- 
fore the daughter looses any appreciable velocity due to interactions with the 
medium (in the superallowed decay of 32 Ar, the daughter travels < 2 x 10~ 2 
A before the proton is emitted). 

• because proton decays in these cases violate isospin symmetry the time scale 
for proton decay is long enough so the widths are much narrower than the 
lepton-recoil broadening (in the superallowed decay of 32 Ar, the daughter state 
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has a width of ss 20 eV while the lepton-recoil broadening has a full-width at 
half-maximum of ~ 25 keV). 

2. Limits on Scalar Weak Interactions from 32 Ar (i + Decay 

We performed an experiment at ISOLDE that measured with high precision 
the energies of protons following the + — > + f3 + decay of 32 Ar. Except for this 
experiment and a previous ISOLDE study of 32 Ar by Schardt and Riisager ||, 
there are no other precise determinations of e-v correlations in pure Fermi tran- 
sitions, which explains why, prior to our new result, limits on scalar couplings 
were rather poor in comparison to those on tensor currents Q. We here describe 
the salient features of our experiment and the extracted limits on scalar currents. 
More details on the 32 Ar experiment can be found in Ref. 

Critical challenges for this kind of experiment are: 

1. obtaining an intense and pure source of radioactivity, 

2. eliminating proton-/? + summing, which distorts the shape of the proton peak, 

3. and optimizing the energy resolution of the proton counter. 

ISOLDE solved problem 1 by producing very pure beams of both 32 Ar and 33 Ar 
from a CaO target and a plasma ion source, providing an average of ~ 94 32 Ar's/s 
and « 3900 33 Ar's/s on our catcher foil over the 9-day-long run. We solved 
problem 2 by immersing our detection system in a 3.5 T magnetic field. The 
32 Ar beam from ISOLDE was stopped in a ~23 /xg/cm 2 C foil at 45 degrees to 
the beam. Our proton detectors were located at ±90 degrees with respect to the 
beam and at about 1.6 cm from the beam spot. In the 3.5 T field superallowed 
/9's had a maximum radius of ~ 0.55 cm, while protons had ~ 7.14 cm. We 
solved problem 3 by using cooled PIN-diode proton detectors (« 0.9 x 0.9 cm) 
and temperature-controlled electronics to obtain a proton energy resolution of 
~ 4.5 keV (~ 3.0 keV electronic noise). This greatly enhanced our sensitivity to 
the e-v correlation. 

Fig. [I] shows Monte-Carlo predictions for the shape of the proton peak as- 
suming the proton detector had infinitely good energy resolution. If one assumes 
b = in Eq. | one can extract a by producing linear combinations of these two 
shapes, folding them with the detector response function (assumed to be two 
exponentials with adjustable tails and areas, convoluted with a Gaussian of ad- 
justable width) and determining the values of the parameters that minimize \ 2 • 
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Figure 1. Intrinsic shapes of the delayed proton group from 32 Ar + — > + decay for a = +1, b = 
(unshaded curve) and a = — 1, 6 = (light-shaded curve). 

Fig. |2] shows our best fit to a subset of our data containing approximately 1/10 
of the statistics. The very good energy resolution obtained in this study also 
allowed us to extract valuable spectroscopic information on 32 CI. The fit was 
performed using an R-matrix parametrization of the resonances. Reduced total 
widths, T p , and ratios T p i/T p o can be extracted for many of the resonances. For 
the general case when b ^ one has to fold an additional distribution (taking into 
account the m/E term in Eq. |2|). We produced a grid in the Cs, C' s space and, 
for each point, minimized x 2 with respect to the response function parameters. 
The resulting confidence regions are shown in Fig. |3|. We found that our Cs, C' s 
constraints are well reproduced by the single parameter 

a = a/(l + 0.19136) 

(7) 

where a and b are given in Eqs. ||, ||. In other words, replacing m/E by an 
appropriate average reproduces the regions of interest and allows us to quote our 
results in terms of o. Our experiment yields the constraint 

a = 0.9989 ± 0.0052(stat.) ± 0.0036(syst.) 68% c.l. (8) 
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Figure 2. R-matrix fit to the 32 Ar delayed proton spectrum. This spectrum, contains roughly 

1/10 of our data. 



The systematic error was evaluated by redoing the whole analysis under different 
conditions. We found da/dA = —1.2 x 1CP 3 keV -1 where A is the /3-decay 
endpoint; and da/dQ p = —0.9 x 1CT 3 keV -1 , where Q p is the energy of the 
emitted proton. We measured Q p with an uncertainty, 5Q P = ±1.2 keV, by 
alternating between ~ 2 h 32 Ar runs with 10-15 min 33 Ar runs that gave us a 
continuous calibration of the energy scale. The mass of 32 Ar has been determined 
only to within 50 keV ||, which would impose a systematic error of ~ 6% on 
our measurement. Fortunately, as shown in Table p], the masses of all other 
members of the T = 2 isospin multiplet are known with high precision. We use 
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Figure 3. 95% conf. limits on Cs and C's, including statistical and systematic errors. Left 
panel: time-reversal-even couplings. The annulus is from this work, the diagonal band is the 
Fierz interference result of Ref. Q. Right panel: time-reversal-odd couplings. The circles are 
from this work and correspond to phases of Cs and C' s of ±90°, +45° and —45°. The shaded 
oval is the constraint with no assumptions about this phase. The diagonal band is from the 

^-coefficient in 19 Ne decay M. 



Table 1 

Comparison of the measured mass excesses of the lowest T — 2 quintet in A = 32 to predictions 



of the Isospin-Multiplet Mass Equation [P(x 2 



0.73]. 





isobar 


T 3 


M cxp (keV) a 


Mjmme (keV) 


32 Si 


+2 


-24080.9 ± 2.2 


-24081.9 ± 1.4 


32p 


+1 


-19232.88 ± 0.20 b 


-19232.9 ±0.2 


32g 





-13970.98 ±0.41 c 


-13971.1 ±0.4 


32 C1 


-1 


-8296.9 ± 1.2 d 


-8296.6 ± 1.1 


32 Ar 


-2 


-2180 ±50 


-2209.3 ±3.2 



"unless noted otherwise, ground state masses are from Ref. [g]. 
b E x = 5072.44 ± 0.06 keV from Ref. J]. 
C E X = 12045.0 ± 0.4 keV from Ref. |k,0. 

d from delayed proton energy measured here and masses of Ref. 



the Isospin-Multiplet Mass Equation |T(J, M(Ts) = a + W3 + cT| , to obtain 



A = 6087.3 ± 2.2 keV. 



(9) 



As shown in Fig. we obtain an excellent fit to the IMME with P(x 2 : u) 
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Figure 4. Plot of residuals of IMME fit to the T = 2 quintuplet in yl = 32. Top panel: mass 
excesses in MeV. The points show the measured values and the continuous line shows the IMME 

fit. Bottom panel: fit residuals in keV. 



0.71. On the other hand, modifying the IMME by adding a dT% term we obtain 
A = 6086.7 ± 4.9 keV and d = 0.25 ± 0.47 keV, but with a lower probability, 
P(x 2 i v ) = 0.52, indicating that there is no empirical basis for adding a term to 
the IMME. We expect the IMME to work very well for the A = 32 multiplet where 
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the states are relatively well-bound and the Coulomb barriers relatively high. 
Small non-zero d terms (never more significant than 3a) have been observed in 
light nuclei where the T3 = T members are much closer to being unbound (or even 
unbound) and the Coulomb barriers much lower. We here assume 5 A = ±2.2 keV 
which combined with the uncertainty in Q p yields a kinematic systematic error 
5a = ±0.0032. Note that even when d is allowed to vary freely the uncertainty is 
only about twice this value. We also checked the dependence of a on the fitting 
regions of the proton spectra; a 28% variation in the width of the region changed 
a by less than ±0.00055. We examined the dependence of our results on the 
form of the detector response function by re-analysing the data with a single-tail 
response function; by re-analyzing the data assuming that a weak Gamow- Teller 
peak lay under the tail of the 32 Ar superallowed peak; and by simultaneously 
fitting the 33 Ar and 32 Ar superallowed peaks using a common response function. 
From these tests we inferred a line-shape systematic error of 5a = ±0.0016. 

Figure [| compares our results to previous constraints on scalar couplings. 
For scalar interactions with Cs = — C' s so that 6 = 0, our data yield the la 
constraint \Cs\ 2 < 3.6xl0~ 3 . The corresponding lower limit on the mass of scalar 
particles with gauge coupling strength is Ms = \Cs\~ 1 ^ 2 Mw > 4.1Mw- We note 
that data from neutron (3 decay by itself does not place stringent constraints on 
scalar couplings because the measurement on the correlation coefficient is not 
very accurate (5a /a 5%) and neutron decay is sensitive to both scalar and 
tensor couplings. Moreover, because of the larger value of < m/E > (i.e. the 
lower endpoint energy), the circle generated by the equivalent of Eq. ^ has a 
large radius which weakens the neutron constraints. Even when supplemented 
by other data on the GT Fierz interference etc., the neutron constraints are not as 
tight as those from the present work and limits on Fierz interference in + — * + 
transitions. 



3. Isospin mixing Corrections to + — > + Tt values and \\f u< ^\ 

Taken at face value, the V u d matrix element extracted from the Tt values 
of nine + — > + /3-decay transitions implies non-unitarity of the Kobayashi- 
Maskawa matrix [17]. Because of the importance and unexpected nature of this 
conclusion, it is worth reexamining whether any systematic effect could affect the 
Tt values. One possibility concerns the corrections for isospin-symmetry violation 
in the parent and daughter nuclear wave functions. These corrections are usually 
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Figure 5. Comparison of our constraints on scalar couplings with previous work. The light 
shaded area represents constraints from neutron (5 decay alone: i.e. from measurements of a, 
A, B, and t 1 / 2 (12) ■ A slightly darker area shows how these constraints improve when combined 
with measurements of the polarization of /3's from 14 O and 10 C [Q. The darker shaded area 
shows the result of adding to the previous the constraints on Fierz terms from 22 Na |y| and 
the measurement of a in 6 He |l5| . The darkest shaded area shows constraints from our results. 
The narrow area looking like a line at —45° is from constraints on Fierz terms from + — > + 

transitions. 

separated into 'configuration mixing' and 'nucleon overlap' parts; the latter being 
dominant (~ four times the former). Although several authors have 
performed independent calculations that agree reasonably well it would be valu- 
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able to check these calculations on additional transitions in neighboring nuclei. 

The cases we present below provide a good opportunity to check the cal- 
culated corrections; the 'nucleon overlap' corrections are enhanced over those in 
the nine standard cases because the nuclei lie farther from the valley of stability. 

3.1. Isospin Mixing in the Fermi decay of 32 Ar 

The isospin-mixing correction in 32 Ar was calculated by B.A. Brown 21] 



using the SKX Skyrme interaction [E2]. This calculation yields 2.0 ±0.4% where 
the uncertainty is based on previous comparisons of similar calculations to mea- 
surements. The large size of the correction is due to the looser binding of the 
d| and proton states compared to the neutron states. For comparison, the 
average correction for the nine standard cases is ~ 0.41%, while the correction 
for the neighboring decay of 34 C1 is ~ 0.61%. 

In 32 Ar decay there is no mixing with GT transitions so that one can check 
the isospin-mixing correction simply by determining the half- life of the parent, 
and the branching ratio and endpoint of the superallowed transition. One can 
thus extract the J-t value and compare with the prediction. Furthermore, if 
some Fermi strength is diverted into narrow J T = + ; T = 1 levels in the 32 C1 
daughter, it may be possible to identify these transitions because they will have 
a = +1 instead of a = —1/3. 

We expect to determine the 32 Ar half life to ~ 0.2% from our ISOLDE 
experiment (our preliminary value is t]/ 2 = 100.74 ± 0.18 ms) and we hope to 
determine the branching ratio with high precision at an upcoming experiment at 
MSU ||]. 

The extraction of the 32 Ar mass from the IMME presented in section ||| yields 
the endpoint to ~ 2.2 keV, which allows to calculate the phase space factor to 
sa 0.3%. The 32 Ar mass could also be measured at ISOLTRAP pi but due 



to the short 32 Ar half-life it is not clear that one could get enough intensity to 
determine the endpoint to within a few keV. All of the above indicates that one 
could extract the J-t value to ~ 0.41%. 

3.2. Isospin Mixing in the Fermi decay of 3S Ar 

The theoretical isospin-mixing calculations can also be tested by studying 
the superallowed decays of the A = An + 1 nuclei. These are mixed Fermi/GT 
transitions so that one needs to determine the B(F)/B(GT) ratios as well as the 
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Table 2 

Comparison of the measured mass excesses of the lowest T — 3/2 quintet in A — 33 to predictions 
of the Isospin-Multiplet Mass Equation [P(x 2 ,v) = 0.51]. 



isobar 


T 3 


M exp (keV) a 


Mimme (keV) 


33 p 


+3/2 


-26337.7 ± 1.1 


-26337.7+ 1.1 


33 g 


+ 1/2 


-21106.14 + 0.41 6 


-21106.15 ± 0.41 


33 CI 


-1/2 


-15460.2 ± 1.0 C 


-15460.10 ± 1.0 


33 Ar 


-3/2 


-9380 ± 30 


-9399.54 ± 3.4 



"unless noted otherwise, ground state masses are from Ref. 
'excitation energy from Ref. Q and masses of Ref. 0. 
c from 32 S(p,p) resonance energy || and masses of Ref. Q|. 



half-lives, branching ratios and energy releases of the superallowed transitions. 
The B(F)/B(GT) ratio can be obtained from the positron-neutrino correlation. 
As pointed out in Section || the e-v correlation from mixed Fermi-GT transitions 
is not as useful for extracting information on scalar or tensor currents because 
one needs to know the relative amounts of each component. However, one can 
take the existing limits on scalar and tensor currents from other experiments and 
use the e-v correlation to extract the B(GT)/B(F) ratio. For simplicity, we here 
assume the scalar and tensor couplings to be zero. 

We will soon gather the necessary information for the case of 33 Ar. We will 
obtain the half-life of 33 Ar and the absolute branching ratio of the superallowed 



transition from our experiment at MSU [23]. The endpoint can be extracted from 



the Isospin-Multiplet Mass Equation. Table ^ shows the mass excesses for the 
T = 3/2 quartet and the corresponding values of the IMME fit. Here we get: 

A( 33 Ar) = 6060.6 ± 2.6 keV. (10) 

In addition there are plans to determine the 33 Ar mass to ~ 3 keV using 
ISOLTRAP p^l (the mass of the daughter level is known to ~ 1 keV). 

We obtain a from the 33 Ar delayed proton spectrum, a typical example of 
which is shown in Fig ||. Our fit to the ISOLDE data yields 

a ( 33 Ar) = 0.944 ± 0.002(stat.) ± 0.003(syst.), (11) 



which implies B(GT)/B(F) = 0.044 ± 0.002 Our result is in reasonably 

good agreement with the the shell-model calculation of [2J| which predicts 
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Figure 6. R-matrix fit to the 33 Ar delayed proton spectrum. This spectrum, contains roughly 

1/10 of our data. 

B(GT)/B(F) = 0.055, but in strong disagreement with the previous determi- 
nation of Schardt and Riisager who obtained || 

a( 33 Ar) > 1.02 ± 0.04 (12) 

(2a error bars), which can be translated into an upper limit B(GT)/B(F) < 0.015 
which is ~ 4 of their <r's from our value. 

We evaluated the systematic errors in a following the same procedure we 
used for 32 Ar. Using the endpoint deduced from the IMME and the masses of 
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Table |^, combined with the known Q p = 2276.5 db 1.0 keV, and the derivatives 
da/dA = -9.1 x 10~ 4 keV" 1 ; and da/dQ p = -8.5 x 1(T 4 keV" 1 , extracted by re- 
doing the analysis with different values of A and Q p , yields 5a = 0.0023. Varying 
the width of the fitting region by 28% yields variations of 5d ~ 0.0004. Simulta- 
neously fitting the 32 Ar and 33 Ar data with a common detector response function 
yields 5a ~ 0.002. These three uncertainties were combined in quadratures to 
give the total systematic error shown in Eq. 
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The discrepancy between our result in Eq. 11 and that of Schardt and Ri- 
isager in Eq. [l^ is due primarily to differences in the analysis rather than dis- 
agreement of the data itself. Dieter Schardt kindly made the raw data of Ref. || 
available to us and our analysis of their data gave a result essentially consistent 



with Eq. 11 



We now can show the potential value of measuring the J-t value for the 
superallowed transition by imagining that the total strength B(F) + B(GT) has 
been determined to ~ 0.3%. Combining this with our positron-neutrino correla- 
tion measurement would yield B(F) to ~ 0.5%. The predicted [21] isospin-mixing 
correction for 33 Ar is ~ 1.2%. For comparison, the correction in 34 C1 is ~ 0.6%. 
The fact that the correction is enhanced, as in 32 Ar, makes it easier to measure. 
So these measurements could check whether the calculated corrections are accu- 
rate to within ~50%. It should be noted that larger discrepancies (corrections 
should be 0.7% as opposed to « 0.4% |2^]) would be needed to explain away 
the apparent non-unitarity of the Kobayashi-Maskawa matrix. 
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